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This work is motivated by a particular 
scheduling problem that is faced by logis- 
tics centers that perform aircraft mainte- 
nance and modification. Here we concen- 
trate on a single facility (hangar) which is 
equipped with several work stations 
(bays). Specifically, a number of jobs have 
already been scheduled for processing at 
the facility; the starting times, durations, 
and work station assignments for these 
jobs are assumed to be known. We are 
interested in bow best to schedule a num- 
ber of new jobs that the facility will be 
processing in the near future. We first 
develop a mixed integer quadratic pro- 
gramming model (MIQP) for this problem. 



Since the exact solution of this MIQP for- 
mulation is time consuming, we develop a 
heuristic procedure, based on existing bin 
packing techniques. This heuristic is fur- 
ther enhanced by application of certain 
local optrmality conditions. 
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1. Introduction 

The present study was motivated by a problem 
encountered by an aircraft logistics center that houses a 
number of facilities (hangars), each of which contains 
multiple work stations (bays). The center offers its cus- 
tomers a wide variety of aircraft maintenance services 
(jobs), including structural repairs, passenger/freight 
conversions, and engine changes. In this paper we con- 
centrate on a single hangar which is equipped with sev- 
eral bays. 

It is assumed that the company has contractual obli- 
gations for a number of existing jobs, which have to be 
completed in the hangar These existing jobs are char- 
acterized by their specified starting times, bay assign- 
ments, and job durations {spans). Thus, a given bay will 
be idle during the time periods between existing jobs. 
We are interested in how best to schedule a number of 
new jobs into these vacant spaces. Specifically, the 



company wishes to construct a manufacturing schedule 
that can be flexible in accommodating future incoming 
jobs. 

First we present a mixed integer quadratic program- 
ming (MIQP) approach to model this scheduling prob- 
lem. Solution of this MIQP gives a schedule that opti- 
mizes the overall facility utilization, while providing 
increased flexibility to accommodate future jobs. Next 
we develop local optimality conditions for the common 
types of changes that can occur to a particular schedule 
(job switches). These conditions enable us to better 
understand optimal solutions and also to conduct sensi- 
tivity analyses. 

The exact MIQP formulation of the scheduling prob- 
lem has mn binary variables, where m denotes the num- 
ber of vacant spaces and n denotes the number of new 
jobs. Since there are 2""" possible choices for the binary 
variables, exact solution of our MIQP model is possible 
only for very small problems. Consequently we devel- 



103 



Volume 111, Number 2, March- April 2006 

Journal of Research of the National Institute of Standards and Technology 



op a heuristic procedure by considering the vacant 
spaces as bins with varying capacities and applying 
some modified bin packing techniques. We can also 
incorporate the local optimality conditions to improve 
this basic heuristic approach. 



2. Optimization Model 

First we consider a very simplified scenario in order 
to motivate our optimization model. Suppose there is 
only one bay and, relative to the jobs already scheduled 
for that bay, there are three vacant spaces (bins) B^, B2, 
S3 of lengths (capacities) 8 d, 12 d, 15 d. Also suppose 
that three new jobs Jy J2, J3 with spans of 6 d, 7 d, 10 d 
are to be assigned to the bay. Several ways of assigning 
the three new jobs to the three bins are depicted in Fig. 
1 . Intuitively, it seems that assignment A^ is the most 
desirable of the three, because it provides the largest 
contiguous space remaining in a bin (9 d) after assign- 
ing the three jobs. For example, by adopting assign- 
ment A}, the center is capable of handling a future job 
that requires a span of 9 d or handling two additional 
jobs with spans of 6 d and 3 d, whereas the other two 
assignments cannot. Here we observe that creating a 
large amount of space in a bin (or in a few bins) after 
assigning the jobs is generally preferred over small 
amounts of space that are scattered over many bins. 

To be more precise, let us define the residual capac- 
ity of a bin as the capacity that remains after assigning 
a job (or jobs) to that particular bin. Since the sum of 
the residuals is fixed (in this example totaling 12) for 
all assignments, we consider the sum of squared resid- 
uals for these assignments. Assignment Ai has residual 
capacities 8, 2, 2 with sum of squares 64 + 4 + 4 = 72; 
similarly, A2 has residual capacities 2, 5, 5 with sum of 
squares 54, while A^ has residual capacities 1 , 2, 9 with 
sum of squares 86. We observe that assignment A^ has 
the largest sum of squares, as a result of the contribu- 



«, 



li^ 



tion from the highest residual capacity term. In general, 
we argue that an assignment with the largest sum of 
squares should be preferred because it leaves a higher 
residual capacity in a bin (or bins) instead of smaller 
residual capacities that are scattered over several bins. 
In turn this gives added flexibility to the company in 
accommodating future incoming jobs. We are thus 
motivated to use the maximization of the sum of 
squared residuals as our objective criterion. Since the 
sum of the residuals is fixed, we are equivalently max- 
imizing the variance of the residuals. 

2.1 Formulation 

To formulate this problem, suppose that n new jobs 
are to be scheduled using m existing bins, with Ci, being 
the capacity of bin B,^. Also suppose that job J, has span 
Sj. Define the binary variable x^ = 1 if job J/ is assigned 
to bin Sj, and Xn, = otherwise. Also let z^ denote the 
residual capacity of bin B,^. Then our optimization 
model for the scheduling problem has the form: 



max z = y^Zj^ 

k=l 
n 

Yi^^k+h=Ct, k^\,...,m 

m 
k=l 

x,,e{0,l}, z, >0. 

The above optimization model has m + n equality 
constraints involving mn binary variables and m contin- 
uous variables. The first m equality constraints define 
the residual capacity variables z^; the nonnegativity of 
zi^ ensures that the set of jobs assigned to bin B,^ should 
not exceed that bin's capacity c^. The remaining n 
equality constraints require that each job should be 
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Fig. 1. Assignments ^ J, A2, A^. 



104 



Volume HI, Number 2, March- April 2006 

Journal of Research of the National Institute of Standards and Technology 



assigned to exactly one bin. Since we are optimizing a 
quadratic objective function over the (linear) constraint 
set, this model can be classified as a mixed integer 
quadratic programming (MIQP) problem. Since the 
objective function is a sum of squares, this problem is 
a type of convex maximization problem, a very difficult 
type of mixed integer optimization problem. 

Thus far we have not considered the starting time i^ 
associated with each bin. We can easily incorporate 
such starting times into our MIQP model by simply 
modifying the objective function: 



maxz(J3) = £(z,+i34)^ 



(1) 



where j3 is a nonnegative valued parameter When p = 
0, we recover our original model in which maximizing 
the sum of squared residuals ("good fit") is the sole cri- 
terion. However, by increasing the value of j3, we 
emphasize the contribution of the starting time and thus 
encourage new jobs to be assigned earlier in time. The 
parameter p can be interpreted as the increased impor- 
tance that the decision maker would like to give to ear- 
lier completion of jobs. It is not hard to verify that z{p) 
is a piecewise quadratic, increasing convex function of 

A 

The above model can be easily applied to situations 
in which there are multiple bays within a hangar 
Indeed for the purposes of the optimization model 
MIQP, it is only necessary to maintain a list of bins with 
their capacities and starting times, without regard to the 
particular bay associated with each bin. In the multiple 
bay context, the parameter j3 can still be interpreted as 
the relative importance of a good fit versus an assign- 
ment of jobs to bins occurring earlier in time. 

2.2 Examples 

In order to get some insights into the proposed MIQP 
model we present two numerical examples. We have 
used the global optimization package LINGO [4] to 



solve these MIQP problems exactly. Example 1 
involves m = 5 bins and two bays: the first three bins B^, 
B2, Bj are associated with bay 1 and the other two bins 
S4, B^ with bay 2. Capacities and starting times for 
these bins are given in Table 1 . Four new jobs Ji, J2, J3, 
J4 are to be assigned; these jobs have the span times 5, 
given in Table 2. Using LINGO we obtained the exact 
solution of the MIQP model for various values of p. It 
turned out that there were only three sets of optimal 
solutions as the parameter P was varied. For < P < 
. 1 1 , the optimal assignment ^ 1 allocates j obs Ji , J2, J3, 
J4 to bins B2, B^, B^, B^ respectively. A schematic illus- 
tration of assignment^! is given in Fig. 2. 

Table 1. Bin capacities and starting times for Example 1 
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Table 2. 


Span times of new jobs for Example 1 
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For the range 0.11 <P< 0.4, the optimal assignment 
A2 allocates jobs Ji, J2, Jy, Jn to bins S4, B2, B^, B^ 
respectively. Thus the only modification that occurs to 
the optimal job schedule, in changing from assignment 
A^ to assignment A2, is that jobs Jj, Jj interchange posi- 
tions between bins B2, B^. Over the final range 0.4 < P 
< 00, the optimal assignment A^ allocates jobs Jj, J2, J3, 
J4 to bins S4, B2, B^, B^ respectively. The only difference 
between assignment A2 and assignment A-^ is the move- 
ment of job J3 from bin S3 to bin B^. This behavior is 
consistent with our earlier observation that increasing P 
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Fig. 2. Assigimient^i, optimal for < j8< 0.11. 
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favors moving jobs to bins that occur earlier in time: S5 
has a starting time of 28, earher than the starting time 
43 for S3. 

Example 2 involves m = 9 bins and three bays, each 
containing three bins. We are interested in assigning 
five new job Jj, ...., J5 to these bins. The data for this 
problem are given in Tables 3 and 4. Here we also 
investigated the changes in optimal solutions as the 
parameter fi is varied. For this problem, there were 
eight sets of optimal solutions as the parameter fi was 
varied over the ranges: [0, 0.004], [0.004, 0.033], 
[0.033, 0.05], [0.05, 0.079], [0.079, 0.15], [0.15, 1], 
[1, 3], [3, 00]. As one illustration. Fig. 3 presents the 
optimal solution A^ over the first of these ranges. In 
assignment A^, jobs Ji, J4, J5 are assigned to bins S3, B^, 
B^ whereas jobs J2, J3 are assigned to bin Bg. 



Table 3. Bin capacities and 


starting times for Exj 
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Table 4. Span times of new 


jobs for Example 2 
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3. Local Optimality Conditions 

In the previous section, we observed several trans- 
formations that occurred in changing from one optimal 
schedule to another optimal schedule. Specifically, in 
changing from A2 to A^, a job (J3) moves from one bin 
to another bin, a transformation we call a move. In 
changing from Ai to A2, two jobs (Jy Jj) interchange 
positions between two bins, termed a swap. In addition, 
another simple transformation involves transferring 
jobs from several bins to a common bin {grouping) or 
vice versa (ungrouping). Collectively we call these 
types of changes, and combinations thereof, 70Z) switch- 
es. In this section we will first study the individual 
effects of each type of job switch upon a given assign- 
ment. This enables us to develop local optimality con- 
ditions for these types of job switches, relative to a 
given assignment. Using these local optimality condi- 
tions, it is possible to identify job switches that can 
improve a given feasible solution. Moreover, these 
local optimality conditions can aid us in carrying out a 
sensitivity analysis with respect to the parameter p. 

3.1 Move 

Consider the two assignments A and B shown in Fig. 

4. In assignment A, job Jj is currently assigned to bin S„, 
whereas in assignment B, job J, has been moved to bin 
B^. Otherwise the two assignments are identical. The 
starting times of bins S„ and B^ are 4 and t^ respective- 
ly, and the span time of job J, is 5,. Define a„ (resp. a„) 
as the capacity utilized by the other jobs that are 
already assigned to S„ (resp. B^), and r„ (resp. r„) as the 
remaining capacity of bin S„ (resp. B^) before the 
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Fig. 3. Assignment A^, optimal for < j8 < 0.004. 
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Fig. 4. Job J I moves from bin S„ to bin B^. 



assignment of job J,. The residual capacities of bins S„ 
and B, in assignment A are defined as z„ and z„ respec- 
tively, whereas the corresponding residual capacities of 
bins S„ and B^ in assignment B are defined as z'^ and z[. 

1 nUS Zy Ty Sl, Zy ly^ Z^ Ty , Z^ V^ S'^. 

Let Z^, Z5 denote the objective function values in Eq. 
(1) for assignments A,B. Since the only difference 
between the two assignments is the movement of job J-, 
from bin S„ to bin B„ we can ignore the other contribu- 
tions in the objective function values Z^ and Zg, writing 

Z, = (z„ + ^O' + i?. + PKf = (r„ - s, + PQ' + (r, + PQ' 



If we define j'„ = r„ + ^4, y^ = r^ + iit„ then the above 
expressions simplify to 

2^ = (y«-'5,-)'+/v, ZB=yi + {y,-s^^. 

Consequently AZ = Z^-Zg = 2(y„ -jj^, is the change 
in objective function (1). Assignment^ will be locally 
optimal with respect to the move in Fig. 4 if and only if 
AZ > 0; since 5, > this is equivalent to y^- y^> 0. 
However, if % -ja < then it is beneficial to move job 
Ji from bin S„ to bin B„ assuming this move is feasible. 

It will be convenient to use an alternative representa- 
tion to depict the move shown in Fig. 4. Namely, we 
denote the bins S„ and B^ by vertices u and v, and we 
use a directed edge / to denote the movement of job J, 
from bin S„ to bin B^. This is depicted in Fig. 5. 



© — KD 



Fig. 5. Graphical representation of a move. 



To illustrate how this local optimality condition can 
be applied, consider Example 1 and assignment A2, 
which allocates jobs Ji, J-i, J^, Jn to bins B^, B2, B^, B^. 
The only feasible moves involve transferring a job from 
one of the bins B^,..., B^ to bin B^. For example, assign- 
ment A2 remains locally optimal for moving job J^ from 
bin S3 to Bs if and only if < j'j - J3 = (16 + 28)8) - 
(10 + 43)8) = 6 - 15)8, so that p < 0.4. The other condi- 
tions < J5 - ji, < j'5 - j'2> ^ J5 - J4 hold automati- 
cally since P > 0. The result is that assignment A2 is 
locally optimal under moves for all P < 0.4, consistent 
with the results obtained earlier from LINGO. 



Zg = iz: + PQ' + iz: + Pt^y = ir, + PQ' + ir,-s, + PQ' 3.2 Swap 



We now consider the swapping of jobs, relative to 
two assignments A and B. In assignment A, jobs J, and 
Ji, are assigned to bins S„ and B^ respectively, whereas 
in assignment B, jobs Jj and Jj are assigned to bins B^ 
and S„. Therefore the only difference between the two 
assignments is that jobs J, and J^ are swapped between 
bins S„ and B^. This is illustrated in Fig. 6. Again define 
a„ (resp. a„) as the capacity utilized by the other jobs 
that are assigned to S„ (resp. B^), and r„ (resp. r„) as the 
remaining capacity of bin S„ (resp. B^) before the 
assignment of jobs Jj and J,,. The residual capacities of 
bins S„ and B^ in assignment A are defined as z„ and z„ 
respectively, whereas the corresponding residual capac- 
ities of bins S„ and B^ in assignment B are defined as z^ 
and z'. Thus z^ = r^- s^, z, = r,- s^, z^ = r^- s,,, z', = r,- 
s,. 



(eC> 



Fig. 6. A swap of jobs J,- and Jj. 
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As before, let Z^, Z^ be the objective function values 
for assignments A, B. Since the only difference between 
the two assignments is the swapping of jobs J, and J^ 
between bins S„ and B„ we can ignore the other contri- 
butions in the objective function values Z^ and Zg, writ- 
ing 

Z, = (z„ + pQ^ + (z, + PQ' = (r„ - ., + PQ' 
+ (r, -st + pKf = (yu- sd' + (y. - s,f 

Z, = (z'„ + PQ' + (z: + pQ' = (r„ -s, + pt^f 

+ (r, - 5,. + pQ^ = (y„ - s,y + (y, - 5/ 

Therefore AZ = Z^- Zg = 2(y,- jj^,- + 2(y„ - y,)si, = 
2(y„ - j'„)(5, - Si,) and assignment ^ is locally optimal 
with respect to the indicated swap if and only if AZ> 0. 
Otherwise it is beneficial to swap jobs J^ and Jj, assum- 
ing that such a swap is feasible. 

By way of illustration, consider assignment Ai of 
Example 1, shown in Fig. 2. The only feasible swaps 
involve exchanging jobs Ji and J2 or exchanging jobs J3 
and J4. Assignment Ai remains locally optimal for 
swapping job Ji in bin B2 with job J2 in ^4 if and only if 
(y4 -y2)isi - S2) > 0. Simplifying produces < [(14 + 
4j3) - (12 + 22/3)1(12 - 11) = 2 - 18/3 or jS < }. 
Swapping jobs J3 and J4 between bins B^ and Bi impos- 
es the condition (y^ -yi)is4 - ^3) > which simplifies to 
37j3 > 0, which clearly holds for j3 > 0. In summary, 
assignment Ai is locally optimal for swaps over the 
range < j3 < 7, consistent with the results obtained ear- 
lier from LINGO. 

3.3 Grouping and Ungrouping 

Another simple job switch involves the grouping of 
jobs, in which/) jobs are moved from separate bins to a 
new bin. Specifically, suppose that in assignment A, 
jobs Ji, ..., Jp are assigned to bins S„ , ..., S„ , whereas in 
assignment B, they are all grouped together in the sin- 
gle bin B,. This transformation is illustrated in Fig. 7. 
Here we can define the utilized capacity of all bins after 
the assignment of jobs other than J^, ..., J^,, as well as 
the residual capacities of bins S„ , ..., S„ , in a manner 
similar to that done previously. 




Fig. 7. Grouping of jobs Jj, ..., J^. 

The associated objective function values for these 
assignments are given by 

2. =Sk, +ptj+iz^+piy =X(^, -^ +iiif 

ZB=t (< +K fH<+iiKf=ti ^, + Pi^ y 

1=1 1=1 

+ (';-S^,+R)' 

1=1 

Simplification then produces 



AZ = Z.-Z„=2 






An analogous development produces the following 
expression when/? jobs Jy ..., J^, are moved from a sin- 
gle bin S„ top separate bins B^, ...,B^: 



AZ = Z,-Z„ 



X^-^v, "•K,)•^.+X•^.•^y 



Notice the change in sign in the last summation from 
negative to positive for the ungrouping AZ, compared 
with the grouping AZ. 

To illustrate the ungrouping of jobs, consider assign- 
ment A^ in Example 2, shown in Fig. 3. Suppose that 
jobs J2, J3 are moved from bin B^ to bins Bj, B^ respec- 
tively. (This is a feasible ungrouping of jobs.) Then 

AZ = 2{(y-,-yg)s2 + (yi -y9)Si + ^2^3] 

= 2([(15 + 3i3) - (21 + 51i3)]ll + [(17 + 5p) 

-(21 +51j3)]10+ll ■ 10) = 2(4-988j3). 
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Consequently, the proposed ungrouping will be advan- 
tageous when AZ < or j3 > ^. 

3.4 More Complex Job Switches 

This section considers combining the previous trans- 
formations (moves, swaps, groupings, and ungroup- 
ings) into more complex job switches. To motivate the 
general case, we first consider an example involving 
eight jobs and five bins. Namely, the transformations 
that occur, in changing from the current assignment A 
to the new assignment B are as follows: jobs Jj, J5 swap 
between bins B^, B^; jobs J2, J3 ungroup from bin B^ to 
bins Bj, B^\ job J4 moves from Sj to B^ while job J(, 
moves from B^ to S^,; and jobs J-,, Jg ungroup from B^ to 
B^, B^. See Fig. 8 for a graphical representation of this 
more complex rearrangement. Following the previous 
development, we obtain 

ZA = (ya-Si-S2- s^f + (Vi - 54)' + 0, - Ssf +(yd- Sef 

+ (ye-ST-Ssy 
ZB=(ya- Ssf + Oi - ^6 - Sjf + (yc-Sl-S4- Ssf 

+ (yd-s2f + (ye-s3f 

AZ = 2[(y, - yjsi + (yj- j'Js'z + (v. - ya)Si + 0. - >'6>4 

+ iya -yc)Si + (yt - j'A + (y* -ye)si + (yc - j'.^s + siS2 

+ 5'25'3 + 5'j5'3 + S'yS'g — SiSn — 5'45'g — SiS^ — Sf^-i\ . 




Fig. 8. Graphical representation of a more complicated scenario. 

Notice that for each directed edge (/, j) representing 
the movement of job Jf. in Fig. 8 there is a term 2(y, - 
y^s^ in AZ. There are also terms in AZ to represent the 
ungrouping and grouping of jobs. Jobs Jj, J2, J3 
ungroup at vertex a, giving rise to the product terms 
S1S2 + S2S3 + SiS^; likewise the ungrouping of jobs J-,, Jg 
at vertex e gives the term SjS^. On the other hand, jobs 



Ji, J4, Jg group at vertex c, producing the (negative) 
term -(5'i5'4 + s^Ss + SiSg), while jobs J^, Jj group at ver- 
tex b, giving -s^j. 

In general suppose that the directed graph G = {V,E) 
represents the specified changes to a current assign- 
ment, where the set Fof vertices represents the bins and 
the set E of directed edges represents the movement of 
jobs between bins. Denote the span time of the job 
associated with edge (i, J) by Sy. Also let r^(/) denote 
the set of edges leaving vertex / and let r^(/) denote the 
set of edges entering vertex ;. Then the expression for 
AZ becomes 



-Y,Y,^^aSi-a,bsY- (/), a it b). 



4. Heuristic Procedures 

In Sec. 2.1, we modeled the scheduling problem as 
an MIQP with mn binary variables and m continuous 
variables, where m is the number of bins and n is the 
number of new jobs. Since the number of possible 
choices 2""" for the binary variables rapidly becomes 
large, even for small m and n, the exact mathematical 
solution of the MIQP model is very time consuming. 
Therefore we consider heuristic solution approaches in 
this section, rather than exact procedures. 

As noted earlier the vacant time interval between any 
two existing jobs can be considered a bin. Since the 
intervals between existing jobs are of different length, 
the corresponding bins have variable size. The underly- 
ing problem is then to assign the new jobs to these bins 
in such a way that the set of new jobs assigned to a bin 
fits within that bin's capacity; this is to be done in an 
"optimal" fashion. Thus our scheduling problem can be 
viewed as a variable-sized bin packing problem [3,5] 
with an unusual type of objective function. In this sec- 
tion we will first develop a simple heuristic procedure 
based on existing bin packing techniques. Next we will 
apply the local optimality conditions developed in Sec. 
3 to improve this bin packing heuristic. 

4.1 Bin Packing Heuristic 

We begin by reviewing two bin packing algorithms 
that are well known in the literature [2]. They are the 
First Fit (FF) algorithm and the Best Fit (BF) algo- 
rithm. The objective of such standard bin packing algo- 
rithms is to minimize the number of bins that are need- 
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ed to pack a given set of items. The FF algorithm 
assigns the next job into the lowest indexed bin into 
which it will fit. If the next job does not fit into any 
existing bin then we open a new bin and place the next 
job in the new bin. The BF algorithm is similar to the 
FF algorithm, except that it assigns the next job into 
that bin which will leave the smallest residual capacity 
after the assignment. 

In our heuristic algorithm we first order the jobs in 
order of nonincreasing span times. There is an intuitive 
appeal to ordering the span times in this manner: we 
assign the jobs with higher span times first and hope 
that we can accommodate the jobs with smaller span 
times using the spaces that remain. A similar strategy is 
adopted in existing bin packing algorithms [1,6]. 
Before assigning the next job (in order of nonincreas- 
ing span) to a bin, we first identify the bins into which 
the job can fit. Among these candidate bins, we select a 
bin, say bin S„, having the minimum value of z„ + ^4. 
As shown in the proposition below, this procedure 
enables us to (locally) improve the current solution as 
much as possible at the next step. Next we assign the 
new job to the selected bin and update the residual 
capacity of that bin. We follow this procedure until all 
jobs have been assigned. We now summarize more for- 
mally the steps of this heuristic algorithm. 

binjjacking 

1 . Order the jobs by nonincreasing span time. 

2. For the next job J^ in order, with span s^, con- 
sider only those bins with residual capacity at 
least s^. Among such bins, select bin S„ to 
minimize z„ + ^4. 

3. Assign job J^ to bin S„ and update z„<^ z„- 

St 

4. If there are additional jobs to be processed, 
go to Step 2. 

Here ^8 > is the same parameter introduced in Sec. 
2.1. When ^ = we assign the next job J^ into a bin 
(into which it fits) having the minimum current residual 
capacity. Since this will also leave the smallest residual 
capacity after the assignment, our heuristic algorithm is 
similar to the BF algorithm for (5=0. When (5 is large, 
we assign jobs to bins in order of increasing starting 
time. In other words, jobs are assigned to bins that 
occur earliest in time. This is analogous to assigning a 
job to the lowest indexed bin into which it will fit, when 
bins are ordered by starting times. So when ji is large, 
our heuristic algorithm behaves similar to the FF algo- 
rithm. In this way we have blended both BF and FF into 
our particular bin packing heuristic. We now demon- 



strate that our heuristic performs the locally "best" 
assignment for the current job at each step of the algo- 
rithm. 

Proposition. The bin packing heuristic locally 
improves the objective function by as much as possible 
at the next step. 

Proof. Suppose that job Jj, with span s,, is to be 
assigned and that bins B^, ..., Bp are the bins into which 
Jf, can fit. Let z^, ..., z^ be the current residual capacities 
of these bins. Select bin B„ such that 



z„ + ^84 = min{zi + Pt„ ..., Zp + ptj,} . 



(2) 



Let Zj be the objective function value (1) obtained by 
assigning job J^ to bin Bj at the next step. We claim that 
Z„ > Zj for all 1 <j <p. Let a denote the contribution to 
the objective function from all bins other than S„ and Bj 
in the current assignment. Thus 

Z, = a + {z,-s, + Xf + {Zj + pt)^ 

Zj=a+ (z„ + pt:f + {Zj -s, + Ptf 

giving 

Zu - Zj = -2z^i - iptji, + 2zjSt + iptjS^ 

= 2[(z, + A,)-(z„ + A„)te>0, 

where the final inequality follows from Eq. (2). n 

We illustrate the bin packing heuristic using 
Example 1 when fi = 0.3. Ordering the four jobs by 
nonincreasing span times places them in the sequence 
Ji, J2, Jn, Jy To begin, job Jj can fit into bins B2, B^, B^, 
which have residual capacities 12, 14, 16 and starting 
times 22, 4, 28. Since min{12 + 0.3(22), 14 + 0.3(4), 
16 + 0.3(28)} = 15.2 is achieved for bin B^, we assign 
job Ji to bin B^ and update the residual capacity of bin 
B4 to Z4 = 2. 

We next select (in order) job J2, which can fit into 
bins B2, B^, with residual capacities 12, 16. Since 
min{12 + 0.3(22), 16 + 0.3(28)} = 18.6 is achieved for 
bin B2, we assign job Jj to bin B2 and update the resid- 
ual capacity of bin Sj to Zj = 1. Continuing in this fash- 
ion, job J4 is assigned to bin Bi and job J3 is assigned to 
bin By Thus the heuristic assigns jobs Jy J2, Jy J^ to 
bins S4, B2, By Bi respectively. This assignment is 
identical to the optimal assignment A2 found for the 
range 0.11 < ^8 < 0.4. In fact for Example 1 the heuris- 
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tic produces optimal assignments over all three ranges 
for the parameter /?. 

For Example 2 the MIQP model gives eight optimal 
assignments corresponding to eight ranges of the 
parameter /?. The assignments obtained by the heuristic 
procedure are identical to the optimal assignments 
obtained by the MIQP model in seven of these eight 
ranges. The only difference occurs over the range 
< P < 0.004 for which the optimal assignment A^ is 
shown in Fig. 3; it assigns jobs Jy ..., J^ to bins B^, B^, 
Bg, B^, B^. By contrast, the heuristic procedure obtains 
a different assignment A'^, in which Jy ..., J^ are 
assigned to B^, Bj, By B^, B^. If however we perform in 
^4 a grouping of jobs J2, J3 from bins Bj, B^ to bin B^, 
then the change in objective function value is 

AZ = 2[(y, -yT)s2 + (y.-y^s, - s^s,] = 2(-4 + 988/?). 

For 0<P< 0.004, the term in parentheses above is neg- 
ative so that it is advantageous to perform this group- 
ing. In other words, first applying the heuristic to obtain 
^4 and then using an improving step (grouping) does 
indeed yield the optimal assignment A^ over the range. 
To summarize, for Example 1 the heuristic procedure 
obtained the optimal solution for all ranges. In Example 
2, there was one instance in which the heuristic gave a 
suboptimal solution. However in this case, a single job 
switch (grouping) was sufficient to produce the optimal 
solution. This encouraging success suggests a hybrid 
heuristic procedure that first carries out the bin packing 
algorithm, followed by local improvements using 
selected job switches. In particular, it is straightforward 
to check for improving moves and swaps; it is a bit 
more tedious to evaluate all groupings and ungroupings 
relative to a given assignment. 



5. Conclusions 

In this paper we have considered the scheduling of 
different types of aircraft maintenance programs. In this 
initial study, we concentrated on the simplified case 
where there is only a single hangar First we developed 
an MIQP to model this scheduling problem. The MIQP 
model incorporates a parameter /? that reflects the rela- 
tive importance of a good fit versus assigning of jobs to 
bins occurring earlier in time. For small test problems, 
it is possible to obtain an exact schedule by solving the 
MIQP model. Results from these test problems indicat- 
ed that there were relatively few optimal schedules over 
the range of all possible values of the specified param- 
eter p. We also developed local optimality conditions 



for certain types of job switches, relative to a given 
assignment. The local optimality conditions enable us 
to improve a given feasible job schedule. Because exact 
solution of the MIQP is limited to fairly small 
instances, we developed a simplified heuristic proce- 
dure based on existing bin packing techniques. An area 
for future research is to combine the bin packing 
heuristic with the intelligent application of the local 
optimality conditions. 
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